Comparative studies of the sensitivities of sparse and full geometries of Total-Body PET scanners built from crystals and plastic scintillators.

BACKGROUND: Alongside the benefits of Total-Body imaging modalities, such as higher sensitivity, single-bed position, low dose imaging, etc., their final construction cost prevents worldwide utilization. The main aim of this study is to present a simulation-based comparison of the sensitivities of existing and currently developed tomographs to introduce a cost-efficient solution for constructing a Total-Body PET scanner based on plastic scintillators. METHODS: For the case of this study, eight tomographs based on the uEXPLORER configuration with different scintillator materials (BGO, LYSO), axial field-of-view (97.4 cm and 194.8 cm), and detector configurations (full and sparse) were simulated. In addition, 8 J-PET scanners with different configurations, such as various axial field-of-view (200 cm and 250 cm), different cross sections of plastic scintillator, and multiple numbers of plastic scintillator layers (2, 3, and 4), based on J-PET technology have been simulated by GATE software. Furthermore, Siemens' Biograph Vision has been simulated to compare the results with standard PET scans. Two types of simulations have been performed. The first one with a centrally located source with a diameter of 1 mm and a length of 250 cm, and the second one with the same source inside a water-filled cylindrical phantom with a diameter of 20 cm and a length of 183 cm. RESULTS: With regards to sensitivity, among all the proposed scanners, the ones constructed with BGO crystals give the best performance ([Formula: see text] 350 cps/kBq at the center). The utilization of sparse geometry or LYSO crystals significantly lowers the achievable sensitivity of such systems. The J-PET design gives a similar sensitivity to the sparse LYSO crystal-based detectors while having full detector coverage over the body. Moreover, it provides uniform sensitivity over the body with additional gain on its sides and provides the possibility for high-quality brain imaging. CONCLUSION: Taking into account not only the sensitivity but also the price of Total-Body PET tomographs, which till now was one of the main obstacles in their widespread clinical availability, the J-PET tomography system based on plastic scintillators could be a cost-efficient alternative for Total-Body PET scanners.

Performance assessment of the 2 γpositronium imaging with the total-body PET scanners.

PURPOSE: In living organisms, the positron-electron annihilation (occurring during the PET imaging) proceeds in about 30% via creation of a metastable ortho-positronium atom. In the tissue, due to the pick-off and conversion processes, over 98% of ortho-positronia annihilate into two 511 keV photons. In this article, we assess the feasibility for reconstruction of the mean ortho-positronium lifetime image based on annihilations into two photons. The main objectives of this work include the (i) estimation of the sensitivity of the total-body PET scanners for the ortho-positronium mean lifetime imaging using 2γ annihilations and (ii) estimation of the spatial and time resolution of the ortho-positronium image as a function of the coincidence resolving time (CRT) of the scanner. METHODS: Simulations are conducted assuming that radiopharmaceutical is labeled with 44Sc isotope emitting one positron and one prompt gamma. The image is reconstructed on the basis of triple coincidence events. The ortho-positronium lifetime spectrum is determined for each voxel of the image. Calculations were performed for cases of total-body detectors build of (i) LYSO scintillators as used in the EXPLORER PET and (ii) plastic scintillators as anticipated for the cost-effective total-body J-PET scanner. To assess the spatial and time resolution, the four cases were considered assuming that CRT is equal to 500 ps, 140 ps, 50 ps, and 10 ps. RESULTS: The estimated total-body PET sensitivity for the registration and selection of image forming triple coincidences (2γ+γprompt) is larger by a factor of 13.5 (for LYSO PET) and by factor of 5.2 (for plastic PET) with respect to the sensitivity for the standard 2γ imaging by LYSO PET scanners with AFOV = 20 cm. The spatial resolution of the ortho-positronium image is comparable with the resolution achievable when using TOF-FBP algorithms already for CRT = 50 ps. For the 20-min scan, the resolution better than 20 ps is expected for the mean ortho-positronium lifetime image determination. CONCLUSIONS: Ortho-positronium mean lifetime imaging based on the annihilations into two photons and prompt gamma is shown to be feasible with the advent of the high sensitivity total-body PET systems and time resolution of the order of tens of picoseconds.

Noise reduction using a Bayesian penalized-likelihood reconstruction algorithm on a time-of-flight PET-CT scanner.

PURPOSE: Q.Clear is a block sequential regularized expectation maximization (BSREM) penalized-likelihood reconstruction algorithm for PET. It tries to improve image quality by controlling noise amplification during image reconstruction. In this study, the noise properties of this BSREM were compared to the ordered-subset expectation maximization (OSEM) algorithm for both phantom and patient data acquired on a state-of-the-art PET/CT. METHODS: The NEMA IQ phantom and a whole-body patient study were acquired on a GE DMI 3-rings system in list mode and different datasets with varying noise levels were generated. Phantom data was evaluated using four different contrast ratios. These were reconstructed using BSREM with different ß-factors of 300-3000 and with a clinical setting used for OSEM including point spread function (PSF) and time-of-flight (TOF) information. Contrast recovery (CR), background noise levels (coefficient of variation, COV), and contrast-to-noise ratio (CNR) were used to determine the performance in the phantom data. Findings based on the phantom data were compared with clinical data. For the patient study, the SUV ratio, metabolic active tumor volumes (MATVs), and the signal-to-noise ratio (SNR) were evaluated using the liver as the background region. RESULTS: Based on the phantom data for the same count statistics, BSREM resulted in higher CR and CNR and lower COV than OSEM. The CR of OSEM matches to the CR of BSREM with ß = 750 at high count statistics for 8:1. A similar trend was observed for the ratios 6:1 and 4:1. A dependence on sphere size, counting statistics, and contrast ratio was confirmed by the CNR of the ratio 2:1. BSREM with ß = 750 for 2.5 and 1.0 min acquisition has comparable COV to the 10 and 5.0 min acquisitions using OSEM. This resulted in a noise reduction by a factor of 2-4 when using BSREM instead of OSEM. For the patient data, a similar trend was observed, and SNR was reduced by at least a factor of 2 while preserving contrast. CONCLUSION: The BSREM reconstruction algorithm allowed a noise reduction without a loss of contrast by a factor of 2-4 compared to OSEM reconstructions for all data evaluated. This reduction can be used to lower the injected dose or shorten the acquisition time.

NEMA NU 2-2007 performance characteristics of GE Signa integrated PET/MR for different PET isotopes.

BACKGROUND: Fully integrated PET/MR systems are being used frequently in clinical research and routine. National Electrical Manufacturers Association (NEMA) characterization of these systems is generally done with 18F which is clinically the most relevant PET isotope. However, other PET isotopes, such as 68Ga and 90Y, are gaining clinical importance as they are of specific interest for oncological applications and for follow-up of 90Y-based radionuclide therapy. These isotopes have a complex decay scheme with a variety of prompt gammas in coincidence. 68Ga and 90Y have higher positron energy and, because of the larger positron range, there may be interference with the magnetic field of the MR compared to 18F. Therefore, it is relevant to determine the performance of PET/MR for these clinically relevant and commercially available isotopes. METHODS: NEMA NU 2-2007 performance measurements were performed for characterizing the spatial resolution, sensitivity, image quality, and the accuracy of attenuation and scatter corrections for 18F, 68Ga, and 90Y. Scatter fraction and noise equivalent count rate (NECR) tests were performed using 18F and 68Ga. All phantom data were acquired on the GE Signa integrated PET/MR system, installed in UZ Leuven, Belgium. RESULTS: 18F, 68Ga, and 90Y NEMA performance tests resulted in substantially different system characteristics. In comparison with 18F, the spatial resolution is about 1 mm larger in the axial direction for 68Ga and no significative effect was found for 90Y. The impact of this lower resolution is also visible in the recovery coefficients of the smallest spheres of 68Ga in image quality measurements, where clearly lower values are obtained. For 90Y, the low number of counts leads to a large variability in the image quality measurements. The primary factor for the sensitivity change is the scale factor related to the positron emission fraction. There is also an impact on the peak NECR, which is lower for 68Ga than for 18F and appears at higher activities. CONCLUSIONS: The system performance of GE Signa integrated PET/MR was substantially different, in terms of NEMA spatial resolution, image quality, and NECR for 68Ga and 90Y compared to 18F. But these differences are compensated by the PET/MR scanner technologies and reconstructions methods.

Outcome in patients perceived as receiving excessive care across different ethical climates: a prospective study in 68 intensive care units in Europe and the USA.

PURPOSE: Whether the quality of the ethical climate in the intensive care unit (ICU) improves the identification of patients receiving excessive care and affects patient outcomes is unknown. METHODS: In this prospective observational study, perceptions of excessive care (PECs) by clinicians working in 68 ICUs in Europe and the USA were collected daily during a 28-day period. The quality of the ethical climate in the ICUs was assessed via a validated questionnaire. We compared the combined endpoint (death, not at home or poor quality of life at 1 year) of patients with PECs and the time from PECs until written treatment-limitation decisions (TLDs) and death across the four climates defined via cluster analysis. RESULTS: Of the 4747 eligible clinicians, 2992 (63%) evaluated the ethical climate in their ICU. Of the 321 and 623 patients not admitted for monitoring only in ICUs with a good (n = 12, 18%) and poor (n = 24, 35%) climate, 36 (11%) and 74 (12%), respectively were identified with PECs by at least two clinicians. Of the 35 and 71 identified patients with an available combined endpoint, 100% (95% CI 90.0-1.00) and 85.9% (75.4-92.0) (P = 0.02) attained that endpoint. The risk of death (HR 1.88, 95% CI 1.20-2.92) or receiving a written TLD (HR 2.32, CI 1.11-4.85) in patients with PECs by at least two clinicians was higher in ICUs with a good climate than in those with a poor one. The differences between ICUs with an average climate, with (n = 12, 18%) or without (n = 20, 29%) nursing involvement at the end of life, and ICUs with a poor climate were less obvious but still in favour of the former. CONCLUSION: Enhancing the quality of the ethical climate in the ICU may improve both the identification of patients receiving excessive care and the decision-making process at the end of life.

Boosting the precision of mediation analyses of randomised experiments through covariate adjustment.

Analyses of randomised experiments frequently include attempts to decompose the intention-to-treat effect into a direct and indirect effect, mediated by given intermediaries, with the aim to shed light onto the treatment mechanism. Methods from causal mediation analysis have facilitated this by allowing for arbitrary models for the outcome and the mediator. They thereby generalise the traditional approach to direct and indirect effects, which is essentially limited to linear models. The default maximum likelihood methods make use of a model for the conditional distribution of the mediator, given treatment and baseline covariates, but are prone to bias when that model is misspecified. In randomised experiments, specification of such model can be easily avoided, but at the expense of a sometimes major efficiency loss when those baseline covariates are predictive of the mediator. In this article, we develop a compromise approach: it makes use of a model for the mediator to optimally extract information from the baseline covariate data but is insulated from the impact of misspecification of that model; it achieves this by exploiting the known randomisation probabilities. Simulation studies and the analysis of a randomised study show major efficiency gains and confirm our theoretical findings that the default methods from causal mediation analysis are sometimes, although not always, reasonably robust to model misspecification. Copyright © 2017 John Wiley & Sons, Ltd.

Recent developments in time-of-flight PET.

While the first time-of-flight (TOF)-positron emission tomography (PET) systems were already built in the early 1980s, limited clinical studies were acquired on these scanners. PET was still a research tool, and the available TOF-PET systems were experimental. Due to a combination of low stopping power and limited spatial resolution (caused by limited light output of the scintillators), these systems could not compete with bismuth germanate (BGO)-based PET scanners. Developments on TOF system were limited for about a decade but started again around 2000. The combination of fast photomultipliers, scintillators with high density, modern electronics, and faster computing power for image reconstruction have made it possible to introduce this principle in clinical TOF-PET systems. This paper reviews recent developments in system design, image reconstruction, corrections, and the potential in new applications for TOF-PET. After explaining the basic principles of time-of-flight, the difficulties in detector technology and electronics to obtain a good and stable timing resolution are shortly explained. The available clinical systems and prototypes under development are described in detail. The development of this type of PET scanner also requires modified image reconstruction with accurate modeling and correction methods. The additional dimension introduced by the time difference motivates a shift from sinogram- to listmode-based reconstruction. This reconstruction is however rather slow and therefore rebinning techniques specific for TOF data have been proposed. The main motivation for TOF-PET remains the large potential for image quality improvement and more accurate quantification for a given number of counts. The gain is related to the ratio of object size and spatial extent of the TOF kernel and is therefore particularly relevant for heavy patients, where image quality degrades significantly due to increased attenuation (low counts) and high scatter fractions. The original calculations for the gain were based on analytical methods. Recent publications for iterative reconstruction have shown that it is difficult to quantify TOF gain into one factor. The gain depends on the measured distribution, the location within the object, and the count rate. In a clinical situation, the gain can be used to either increase the standardized uptake value (SUV) or reduce the image acquisition time or administered dose. The localized nature of the TOF kernel makes it possible to utilize local tomography reconstruction or to separate emission from transmission data. The introduction of TOF also improves the joint estimation of transmission and emission images from emission data only. TOF is also interesting for new applications of PET-like isotopes with low branching ratio for positron fraction. The local nature also reduces the need for fine angular sampling, which makes TOF interesting for limited angle situations like breast PET and online dose imaging in proton or hadron therapy. The aim of this review is to introduce the reader in an educational way into the topic of TOF-PET and to give an overview of the benefits and new opportunities in using this additional information.

Evaluation of resistive-plate-chamber-based TOF-PET applied to in-beam particle therapy monitoring.

Particle therapy is a highly conformal radiotherapy technique which reduces the dose deposited to the surrounding normal tissues. In order to fully exploit its advantages, treatment monitoring is necessary to minimize uncertainties related to the dose delivery. Up to now, the only clinically feasible technique for the monitoring of therapeutic irradiation with particle beams is Positron Emission Tomography (PET). In this work we have compared a Resistive Plate Chamber (RPC)-based PET scanner with a scintillation-crystal-based PET scanner for this application. In general, the main advantages of the RPC-PET system are its excellent timing resolution, low cost, and the possibility of building large area systems. We simulated a partial-ring scanner based on an RPC prototype under construction within the Fondazione per Adroterapia Oncologica (TERA). For comparison with the crystal-based PET scanner we have chosen the geometry of a commercially available PET scanner, the Philips Gemini TF. The coincidence time resolution used in the simulations takes into account the current achievable values as well as expected improvements of both technologies. Several scenarios (including patient data) have been simulated to evaluate the performance of different scanners. Initial results have shown that the low sensitivity of the RPC hampers its application to hadron-beam monitoring, which has an intrinsically low positron yield compared to diagnostic PET. In addition, for in-beam PET there is a further data loss due to the partial ring configuration. In order to improve the performance of the RPC-based scanner, an improved version of the RPC detector (modifying the thickness of the gas and glass layers), providing a larger sensitivity, has been simulated and compared with an axially extended version of the crystal-based device. The improved version of the RPC shows better performance than the prototype, but the extended version of the crystal-based PET outperforms all other options.

Optimized light sharing for high-resolution TOF PET detector based on digital silicon photomultipliers.

The majority of current whole-body PET scanners are based on pixelated scintillator arrays with a transverse pixel size of 4 mm. However, recent studies have shown that decreasing the pixel size to 2 mm can significantly improve image spatial resolution. In this study, the performance of Digital Photon Counter (DPC) from Philips Digital Photon Counting (PDPC) was evaluated to determine their potential for high-resolution whole-body time of flight (TOF) PET scanners. Two detector configurations were evaluated. First, the DPC3200-44-22 DPC array was coupled to a LYSO block of 15 × 15 2 × 2 × 22 mm(3) pixels through a 1 mm thick light guide. Due to light sharing among the dies neighbour logic of the DPC was used. In a second setup the same DPC was coupled directly to a scalable 4 × 4 LYSO matrix of 1.9 × 1.9 × 22 mm(3) crystals with a dedicated reflector arrangement allowing for controlled light sharing patterns inside the matrix. With the first approach an average energy resolution of 14.5% and an average CRT of 376 ps were achieved. For the second configuration an average energy resolution of 11% and an average CRT of 295 ps were achieved. Our studies show that the DPC is a suitable photosensor for a high-resolution TOF-PET detector. The dedicated reflector arrangement allows one to achieve better performances than the light guide approach. The count loss, caused by dark counts, is overcome by fitting the matrix size to the size of DPC single die.

Reliability of lithium dilution cardiac output in anaesthetized sheep.

BACKGROUND: Cardiac output (CO) measurement with lithium dilution (COLD) has not been fully validated in sheep using precise ultrasonic flow probe technology (COUFP). Sheep generate important cardiovascular research models and the use of COLD has become more popular in experimental settings. METHODS: Ultrasonic transit-time perivascular flow probes were surgically implanted on the pulmonary artery of 13 sheep. Paired COLD readings were taken at six time points, before and after implantation of a left ventricular assist device (LVAD) and compared with COUFP recorded just after lithium injection. RESULTS: The mean COLD was 5.7 litre min(-1) (range 3.8-9.6 litre min(-1)) and mean COUFP 5.9 litre min(-1) (range 4.0-9.2 litre min(-1)). The bias (standard deviation) was 0.3 (1.0) litre min(-1) [5.1 (16.9)%] and limits of agreement (LOA) were -1.7 to 2.3 litre min(-1) (-28.8 to 39.0%) with a percentage error (PE) of 34.4%. Data to assess trending [rate (95% confidence intervals)] included a 78 (62-93)% concordance rate in the four-quadrant plot (n=27). In the half moon polar plot (n=19), the mean polar angle was +5°, the radial LOA were -49 to +35° and 68 (47-89)% of data points fell within 22.5° of the mean polar angle. Both tests indicated moderate to poor trending ability. CONCLUSION: COLD is not precise when evaluated against COUFP in sheep based on the statistical criteria set, but the results are comparable with previously published animal studies.

In vitro and in vivo imaging characteristics assessment of polymeric coils compared with standard platinum coils for the treatment of intracranial aneurysms.

BACKGROUND AND PURPOSE: Conventional platinum coils cause imaging artifacts that reduce imaging quality and therefore impair imaging interpretation on intraprocedural or noninvasive follow-up imaging. The purpose of this study was to evaluate imaging characteristics and artifact production of polymeric coils compared with standard platinum coils in vitro and in vivo. MATERIALS AND METHODS: Polymeric coils and standard platinum coils were evaluated in vitro with the use of 2 identical silicon aneurysm models coiled with a packing attenuation of 20% each. DSA, flat panel CT, CT, and MR imaging were performed. In vivo evaluation of imaging characteristics of polymeric coils was performed in experimentally created rabbit carotid bifurcation aneurysms. DSA, CT/CTA, and MR imaging were performed after endovascular treatment of the aneurysms. Images were evaluated regarding visibility of individual coils, coil mass, artifact production, and visibility of residual flow within the aneurysm. RESULTS: Overall, in vitro and in vivo imaging showed relevantly reduced artifact production of polymeric coils in all imaging modalities compared with standard platinum coils. Image quality of CT and MR imaging was improved with the use of polymeric coils, which permitted enhanced depiction of individual coil loops and residual aneurysm lumen as well as the peri-aneurysmal area. Remarkably, CT images demonstrated considerably improved image quality with only minor artifacts compared with standard coils. On DSA, polymeric coils showed transparency and allowed visualization of superimposed vessel structures. CONCLUSIONS: This initial experimental study showed improved imaging quality with the use of polymeric coils compared with standard platinum coils in all imaging modalities. This might be advantageous for improved intraprocedural imaging for the detection of complications and posttreatment noninvasive follow-up imaging.

Energy harvesting from the beating heart by a mass imbalance oscillation generator.

Energy-harvesting devices attract wide interest as power supplies of today's medical implants. Their long lifetime will spare patients from repeated surgical interventions. They also offer the opportunity to further miniaturize existing implants such as pacemakers, defibrillators or recorders of bio signals. A mass imbalance oscillation generator, which consists of a clockwork from a commercially available automatic wrist watch, was used as energy harvesting device to convert the kinetic energy from the cardiac wall motion to electrical energy. An MRI-based motion analysis of the left ventricle revealed basal regions to be energetically most favorable for the rotating unbalance of our harvester. A mathematical model was developed as a tool for optimizing the device's configuration. The model was validated by an in vitro experiment where an arm robot accelerated the harvesting device by reproducing the cardiac motion. Furthermore, in an in vivo experiment, the device was affixed onto a sheep heart for 1 h. The generated power in both experiments-in vitro (30 µW) and in vivo (16.7 µW)-is sufficient to power modern pacemakers.

Fast generation of 4D PET-MR data from real dynamic MR acquisitions.

We have implemented and evaluated a framework for simulating simultaneous dynamic PET-MR data using the anatomic and dynamic information from real MR acquisitions. PET radiotracer distribution is simulated by assigning typical FDG uptake values to segmented MR images with manually inserted additional virtual lesions. PET projection data and images are simulated using analytic forward projections (including attenuation and Poisson statistics) implemented within the image reconstruction package STIR. PET image reconstructions are also performed with STIR. The simulation is validated with numerical simulation based on Monte Carlo (GATE) which uses more accurate physical modelling, but has 150× slower computation time compared to the analytic method for ten respiratory positions and is 7000× slower when performing multiple realizations. Results are validated in terms of region of interest mean values and coefficients of variation for 65 million coincidences including scattered events. Although some discrepancy is observed, agreement between the two different simulation methods is good given the statistical noise in the data. In particular, the percentage difference of the mean values is 3.1% for tissue, 17% for the lungs and 18% for a small lesion. The utility of the procedure is demonstrated by simulating realistic PET-MR datasets from multiple volunteers with different breathing patterns. The usefulness of the toolkit will be shown for performance investigations of the reconstruction, motion correction and attenuation correction algorithms for dynamic PET-MR data.

A simple, economical, and effective portable paediatric mock circulatory system.

Ventricular assist devices (VADs) and total artificial hearts have been in development for the last 50 years. Since their inception, simulators of the circulation with different degrees of complexity have been produced to test these devices in vitro. Currently, a new path has been taken with the extensive efforts to develop paediatric VADs, which require totally different design constraints. This paper presents the manufacturing details of an economical simulator of the systemic paediatric circulation. This simulator allows the insertion of a paediatric VAD, includes a pumping ventricle, and is adjustable within the paediatric range. Rather than focusing on complexity and physiological simulation, this simulator is designed to be simple and practical for rapid device testing. The simulator was instrumented with medical sensors and data were acquired under different conditions with and without the new PediaFlowTM paediatric VAD. The VAD was run at different impeller speeds while simulator settings such as vascular resistance and stroke volume were varied. The hydraulic performance of the VAD under pulsatile conditions could be characterized and the magnetic suspension could be tested via manipulations such as cannula clamping. This compact mock loop has proven to be valuable throughout the PediaFlow development process and has the advantage that it is uncomplicated and can be manufactured cheaply. It can be produced by several research groups and the results of different VADs can then be compared easily.

Accurate Monte Carlo modelling of the back compartments of SPECT cameras.

Today, new single photon emission computed tomography (SPECT) reconstruction techniques rely on accurate Monte Carlo (MC) simulations to optimize reconstructed images. However, existing MC scintillation camera models which usually include an accurate description of the collimator and crystal, lack correct implementation of the gamma camera's back compartments. In the case of dual isotope simultaneous acquisition (DISA), where backscattered photons from the highest energy isotope are detected in the imaging energy window of the second isotope, this approximation may induce simulation errors. Here, we investigate the influence of backscatter compartment modelling on the simulation accuracy of high-energy isotopes. Three models of a scintillation camera were simulated: a simple model (SM), composed only of a collimator and a NaI(Tl) crystal; an intermediate model (IM), adding a simplified description of the backscatter compartments to the previous model and a complete model (CM), accurately simulating the materials and geometries of the camera. The camera models were evaluated with point sources ((67)Ga, (99m)Tc, (111)In, (123)I, (131)I and (18)F) in air without a collimator, in air with a collimator and in water with a collimator. In the latter case, sensitivities and point-spread functions (PSFs) simulated in the photopeak window with the IM and CM are close to the measured values (error below 10.5%). In the backscatter energy window, however, the IM and CM overestimate the FWHM of the detected PSF by 52% and 23%, respectively, while the SM underestimates it by 34%. The backscatter peak fluence is also overestimated by 20% and 10% with the IM and CM, respectively, whereas it is underestimated by 60% with the SM. The results show that an accurate description of the backscatter compartments is required for SPECT simulations of high-energy isotopes (above 300 keV) when the backscatter energy window is of interest.

Compression and reconstruction of sorted PET listmode data.

BACKGROUND: In nuclear medicine data can be stored in histogram or listmode format. The most popular histogram format is the planar projection format. Due to the increase in detector blocks, the improved energy resolution and the trends towards time of flight, dynamic and gated imaging, it can be more appropriate to store the data in listmode format. The size of the storage in this format increases linearly with the number of properties (positions, energy, time info) while the histogram format increases exponentially. However, the datasize of listmode data also increases linearly with the number of coincidences. Due to the high number of counts in 3D PET this will lead to very large datasets. Therefore a good compression algorithm for listmode data is very important. METHODS: A sorting and compression method is proposed to reduce the amount of space needed to store the listmode dataset. One event is represented by one number without any information loss compared to the original listmode file. The next step is to sort all events into an array of increasing numbers. These data are compressed by the gzip routine. One of the advantages of 3D PET listmode reconstructions is that they result in a more uniform resolution across the field of view (FOV), which is not always true for other reconstruction algorithms. This improved resolution is shown for the listmode data of a gamma camera operating in PET mode. RESULTS: First the effect of positional accuracy in the listmode dataset is evaluated by comparing resolution in the reconstructions. It is shown that the highest accuracy is not necessary and a significant reduction in the size of the dataset can be obtained prior to lossless compression. A further reduction can be obtained by using the proposed sorting and compression techniques. It is shown that the storage space decreases linearly with the logarithm of the number of coincidences. The compression obtained by different acquisition matrices was compared. Finally it is shown that the 3D listmode reconstruction of sorted listmode data is faster because of improved cache behaviour. The method can be applied to any kind of listmode data. The compression factors will improve when the ratio of measured events to possible events increases.

Geometrical importance sampling and pulse height tallies in GATE.

Monte Carlo simulations are widely used to study the behavior and detection of gamma photons in medical imaging devices. Such simulations are computationally expensive. This is why geometrical importance sampling, a variance reduction technique, was recently incorporated into the GEANT4 Monte Carlo code. In order to use this technique for single photon emission computed tomography (SPECT) imaging, it needed to be made compatible with pulse height tallies. These tallies correspond to the number of detected pulses in distinct energy bins, covering an energy spectrum relevant to SPECT. Since each pulse is the combination of different detector hits, the tally bin is not known until the end of an event. In an analog simulation (without variance reduction) this poses no problems as each detected hit can be stored and the pulse can be calculated at the end of each event. Geometrical importance sampling combined with Russian Roulette however introduces branches into the particle history, which results in a much more complicated pulse calculation. This work describes how pulse height tallies are adjusted to geometrical importance sampling and Russian Roulette within GATE, a medical imaging and simulation application based on GEANT4. The validation of this technique is done through SPECT simulations comparing the analog result with the new method.

In vitro evaluation of the PUCA II intra-arterial LVAD.

The "pulsatile catheter" (PUCA) pump is a minimally invasive intra-arterial left ventricular assist device intended for acute support of critically ill heart failure patients. To assess the hydrodynamic performance of the PUCA II, driven by an Arrow AutoCat IABP driver, we used a (static) mock circulatory system in which the PUCA II was tested at different loading conditions. The PUCA II was subsequently introduced in a (dynamic) cardiovascular simulator (CVS) to mimic actual in vivo operating conditions, with different heart rates and 2 levels of left ventricular (LV) contractility. Mock circulation data shows that PUCA II pump performance is sensitive to afterload, pump rate and preload. CVS data demonstrate that PUCA II provides effective LV unloading and augments diastolic aortic pressure. The contribution of PUCA II to total flow is inversely related to LV contractility and is higher at high heart rates. We conclude that, with the current IABP driver, the PUCA II is most effective in 1:1 mode in left ventricles with low contractility.

Validation of planar and tomographic radionuclide ventriculography by a dynamic ventricular phantom.

Although there is increasing interest in the automatic processing of tomographic radionuclide ventriculography (TRV) studies, validation is mainly limited to a comparison of TRV results with data from planar radionuclide ventriculography (PRV) or gated perfusion single photon emission computed tomography (SPECT). The aim of this study was to use a dynamic physical cardiac phantom to validate the ejection fraction (EF) and volumes from PRV and TRV studies. A new dynamic left ventricular phantom was constructed and used to obtain 21 acquisitions in the planar and tomographic mode. The directly measured volumes and EFs of the phantom during the acquisitions were considered as the gold standard for comparison with TRV and PRV. EFs were calculated from PRV by background-corrected end-diastolic and end-systolic frames. Volumes and EFs were calculated from TRV by region growing with different lower thresholds to search for the optimal threshold. EF from PRV correlated significantly with the real EF (r=0.94, P=0.00). The optimal threshold value for volume calculation from TRV in 336 cases was 50% (r=0.98, P=0.00) yielding the best slope after linear regression. When considering these calculated end-diastolic and end-systolic volumes, EF correlated well (r=0.99, P=0.00) with the real EF, and this correlation was significantly (P=0.04) higher than that of the EF from PRV. Our experiments prove that EF measured by TRV yields more accurate results compared with PRV in dynamic cardiac phantom studies.